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APPROXIMATE  CALCULATION  OF  LOW  FORWARD  SPEED 
SHIP  RADIATION  AND  DIFFRACTION  WAVE  PATTERNS 


1.  INTRODUCTION 

In  recent  years,  with  the  introduction  of  remote  sensing  systems  such  as  Synthetic  Aperture 
Radar  (SAR).  infrared  (IR)  and  microwave  (MR)  radiometry,  there  has  been  considerable  interest  in 
assessing  the  detectability  of  ship  generated  waves  in  die  presence  of  ambient  ocean  waves.  Thus,  the 
present  author  has  conducted  comparisons  (with  typical  random  ocean  waves)  of  the  Kelvin  waves  due 
to  a  nonoscillating  ship  traveling  at  constant  forward  speed  in  the  spatial  [1]  and  spectral  [2]  domains, 
and  of  the  waves  due  to  a  ship  at  zero  forward  speed  excited  by  the  ocean  waves  in  the  spatial 
domain  [3].  The  ship  waves  have  been  calculated  by  using  far  field  asymptotic  approximations  of  the 
Green’s  functions  for  the  Kelvin  wake  [4]  and  zero  speed  [5]  cases.  In  complementary  studies  using 
more  accurate  near  field  formulations,  the  Kelvin  wake  has  been  calculated  by  using  elementary  infin¬ 
ite  fluid  Rankine  [6]  and  the  more  complex  free  surface  Kelvin  [7]  sources. 

In  all  of  the  above  studies,  the  ship  waves  have  been  calculated  for  the  zero  speed  or  zero  fre¬ 
quency  limits  of  the  more  general  case  of  an  oscillating,  translating  ship.  The  principal  reason  is  due 
to  the  well  known  considerable  increase  in  the  complexity  of  the  Green’s  function  for  the  general  case 
[8,9,10].  To  the  author’s  knowledge,  there  has  been  no  extensive  calculation  of  the  actual  ship  waves 
for  this  general  case.  Using  an  ingenious  approach,  whereby  the  ship  forward  speed  is  considered  to 
be  small  enough  so  that  its  squares  may  be  systematically  neglected,  Grekas  [11,12]  derives  a  low 
speed  Green’s  function  which  is  considerably  simpler  than  the  corresponding  function  for  arbitrary 
speed.  He  goes  on  to  obtain  a  relatively  simple  far  field  asymptotic  expression.  This  approach  is 
used  here,  along  with  the  formulation  given  in  [3],  to  calculate  a  series  of  approximate  spatial  realiza¬ 
tions  of  wave  elevations  and  slopes  generated  by  a  low  speed  ship  excited  by  typical  ambient  ocean 
waves.  As  in  [3],  realizations  are  also  calculated  for  the  random  sea  waves  and  superposed  sea  and 
ship  waves  to  assess  the  detectability  of  the  ship  waves. 

The  report  starts  by  describing  Grekas’s  approach.  In  particular,  the  behavior  of  his  far  field 
expression  as  a  function  of  the  dimensionless  parameter  Uu/g  (where  U  is  the  ship  forward  speed,  o> 
is  the  circular  oscillation  frequency,  and  g  is  the  gravity  constant)  is  discussed  in  detail.  The  pro¬ 
cedure  previously  used  in  [3]  for  generalizing  the  far  field  expression,  which  is  for  a  single  source 
oscillating  at  a  single  frequency,  to  the  case  of  the  multiple  sources  modeling  the  ship  hull  and  the 
multiple  frequencies  and  headings  modeling  the  random  sea,  is  briefly  summarized.  Modifications  of 
this  procedure,  which  is  for  zero  speed,  to  account  for  the  present  case  of  low  forward  speeds,  are 
outlined.  Spatial  realizations  of  ship  waves,  sea  waves,  and  superposed  waves  for  Sts  States  2,  5, 
and  8  are  presented  in  the  form  of  three-dimensional  and  contour  plots.  The  detectability  of  the  ship 
waves  are  discussed  as  a  function  of  sea  state,  the  type  of  realization  (elevation  or  slope),  the  form  of 
presentation  (three-dimensional  or  contour)  and  the  upper  bound  assigned  to  the  parameter  Uu/g. 
The  report  concludes  with  a  summary  of  the  major  findings. 

2.  THEORETICAL  FORMULATION  FOR  SINGLE  SOURCE 


2.1  Low  Speed  Green’s  Function 

For  a  submerged  pulsating  source  in  the  presence  of  a  uniform  current  in  the  x-direction,  shown 
in  Fig.  1.  the  Green's  function  C  must  satisfy  the  following  well  known  free  surface  condition 


d2G  ,  „„  d2G  ,  „2  d2G  dG  n 

— -  +  21  +  U-  — 7  +  g  ~ —  =  0  on  z  =  0. 

dr  d,xdt  d.x2  dz 
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If  the  strength  Q  of  the  pulsating  source  is  taken  to  vary  with  time  as  a  simple  harmonic  function 


A ' 


Q  =  Re(e%"!'“-')  (2) 

where  Re  denotes  taking  the  real  part,  then  G  can  be  expressed  by 

G  =  G[  (x,  y,  z;  x',  y',  z')  e~,U!t .  (3) 

Equation  (1)  can  then  be  rewritten  in  terms  of  spatial  variables  as  follows 

,  3G,  _  32G,  3G, 

-  crG,  -  2£/f«  +  G2 - -  +  g  -r-  =  0ont=0.  (4) 

dx  Qx2  fa 

It  is  well  known  that  the  Green’s  function  satisfying  Eq,  (4)  is  troublesome  to  evaluate  [8].  One  of 
the  principal  difficulties  is  that  Eq.  (4)  gives  a  higher  order  pole  appearing  in  the  denominator  in 
integral  representations  of  G  i . 

To  minimize  these  difficulties,  Grekas  [11,12]  makes  the  crucial  assumption  that  the  velocity  U 
is  sufficiently  small  so  that  U 2  may  be  neglected  in  Eq.  (4),  resulting  in  the  following  simpler  free 
surface  condition 

dG2  3G2 

-  -  2Viu  +  g  — —  =  0  on  z  =  0.  (5) 

ox  oz 

The  resulting  Green’s  function  G2,  which  now  has  only  a  simple  pole,  represents  a  drastic  simplifica¬ 
tion  of  the  more  general  function  Gj.  The  detailed  derivation  given  in  [11],  and  summarized  in  [13], 
will  not  be  repeated  here. 

By  using  the  principle  of  stationary  phase,  Grekas  goes  on  to  obtain  the  following  relatively  sim¬ 
ple  algebraic  expression  for  the  far  field  values  of  G2 


Q*  8  f,ns.  /  2*  j!W*cos(90-7)+*oZ'--ir/4] 

~T  2*3*  W  )  v 


where 


R  =  Vfc  -  x'f  +  O'  -  y’)2 


0  is  the  wave  propagation  direction 


y  =  tan  1  (y  -  y ')/{x  -  x')  is  the  spatial  direction 


fie o)  - 


g  +  2uU  cos  dQ 

ujq  (co  IJ  cos  60)z 


g  +  2  I'u  cos  80 


is  the  wavenumber  in  terms  of  the  intrinsic  circular  frequency 


g"(9  o)  = 


-g  cos  {Q.J  -  y)u2 


ig  +  2  u)  U  cos  80)2 
8q  is  the  root  of  the  6  —  7  relationship  given  by 
*'(0o)  =  £  s'n  (7  ~  0<j)  +  2 osU  sin  7  =  0. 


<7e) 


(7f) 


It  should  oe  pointed  out  that  the  expression  for  g  (Bq)  in  [11]  is  erroneously  raised  to  a  power  of  3 
instead  of  2. 


>ne  check  on  the  accuracy  of  the  above  expression  is  to  consider  the  case  when  V  is  set  equal 
to  zero.  By  noting  that  in  this  case 


G2  for  V  =  0  becomes 
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(9) 

This  is  identical  to  the  far  field  expression  for  the  zero  speed  Green’s  function  given  in  [5,8]. 


One  indication  of  the  simplicity  of  the  low  speed  function  is  that,  similar  to  the  V  -  0  case,  for 
a  given  spatial  direction  7  there  is  one  and  only  one  wave,  whose  direction  of  propagation  d$  is  given 
by  Eq.  (7f).  The  difference  is  that  for  U  =  0,  =  7,  which  leads  to  a  set  of  circular  ring  waves 

which  are  centered  around  the  source.  These  simple  waves  may  be  compared  to  the  more  compli¬ 
cated  Kelvin  wake  which  has  zero  or  two  waves  depending  on  7,  and  the  waves  due  to  G,  which 
have  zero  to  five  waves  depending  on  7  and  G«/g. 


In  the  present  work,  it  is  not  the  Green’s  function  itself,  but  the  elevation  f  and  x  and  y  slopes 
ft  an<i  f¥  which  are  of  interest.  Keeping  only  the  terms  which  are  linear  in  G2  in  the  Bemouilli 
equation,  the  elevation  f  due  to  a  single  pulsating  source  is  given  by 


$•  =  Re 


dG, 


+  V 


dG, 


(10) 


Taking  d/dx  of  G2  given  by  Eq,  (6)  would,  in  general,  involve  the  following  three  factors: 
1  R  »  8  "(9 0).  an£l  exp  [ik0  R  cos  (B0  -  7)].  However,  taking  d/dx  of  the  first  two  factors  gives 
higher  powers  of  R  in  the  denominator  which  may  be  neglected  in  the  far  field.  In  taking  d/dx  of 
the  exponential  term,  u  is  of  interest  to  note  the  following  relationship 


x')  cos  0O  +  (y  -  y')  sin  0O 


The  derivative  3/<k  of  the  exponential  term  is  then  simply  given  by 


8  Ik. -ft  cos  (#„  -0)  _  .,  n  *rft  cos  (ft,  -  ffj 

~~~  e  —  ikq  cos  t?0  e 

ax 

Using  Eqs.  (6),  (10),  and  (12),  the  elevation  f  is  then  given  by 

f  =  Re  —  [(o>  -  Uk0  cos  9q)  G2] . 

8 

The  slopes  and  fv  are  conveniently  obtained  from  f  by  again  using  Eq.  (11) 

fjt  -  T"  f  “  Re  (»0«»V> 

<7X 

fy  =  “  f  *  Re  'i£0  sin  0O  f). 

dy 


(12) 


(13) 


(14a) 

(14b) 


Figures  2a,  b,  and  c  respectively  show  sample  plots  of  f  calculated  by  using  Eq.  (13)  for 
U  =  0,  U  <0  (directed  in  the  -x  direction),  and  U  >  0.  Figure  2a  shows  the  expected  circular 
ring  waves  which  are  symmetrical  in  all  four  quadrants,  while  Figs.  2b  and  2c  show  that  in  the  case 
of  j  U  j  >0,  symmetry  is  preserved  only  in  the  lateral  y  direction.  In  the  x  direction,  the  crests 
tend  to  be  dense  upstream  and  sparse  downstream,  in  accordance  with  the  sketch  given  by  Grekas 
[11,12]. 


2.2  Significance  of  the  Parameter  v  -  Uas/g 


Grekas  points  out  that  his  formulation  is  valid  only  for 


f 

l 


V 


Ua  " 


(15) 


This  can  be  seen  from  Eq.  (7d)  which  shows  that  the  wave  number  k0  —  oo(X  —  0)  as  v  — 
well  as  from  Eqs.  (7e)  and  (6)  which  show  that  g"  —  0  (and  hence  G2  -*  oo)  as  v  —  y. 


1/2,  as 
Figures 


3a  to  3e  give  three-dimensional  plots  which  show  the  manner  in  which  the  maximum  wave  amplitudes 
increase  and  the  upstream  wavelengths  decrease  as  v  increases  from  0  to  0.49.  Figure  3e  shows 
that  for  v  ~  0.49,  the  chosen  grid  is  not  fme  enough  to  model  the  short  wavelengths  and  the 
upstream  region  is  plotted  as  a  plateau.  Table  1  shows  the  variation  of  the  ratio  /fwo,  where 
is  the  value  of  for  the  U  =0  case,  for  values  of  v  between  0  and  0.49. 


The  above  shows  that  it  is  important  to  place  an  upper  limit  on  v,  vy,  to  prevent  unrealistically 
large  values  of  f.  (For  example,  it  is  shown  later  that  for  typical  sea  states,  it  is  always  possible  to 
make  the  ship  generated  waves  dominate  the  ocean  waves  by  choosing  vv  sufficiently  close  to  the 
singular  limit  of  0.5.)  Basically,  it  is  of  interest  to  choose  a  value  of  vy  which  is  large  enough  to 
give  an  indication  of  the  effect  of  (7,  without  giving  excessively  large  elevations.  It  was  decided  to 
make  most  of  the  calculations  of  the  ship  waves  using  a  value  of  vv  =  0.25.  It  may  be  noted  that 
this  '•  due  lies  halfway  between  the  possible  limits  of  0  and  0.5,  and  that  at  this  value  the  maximum 
elevation  is  approximately  three  times  that  of  the  U  =  0  case  (see  Table  1),  which  represents  a 
moderate  but  appreciable  modification.  Also,  it  is  well  known  that  for  the  case  of  the  general 
Green's  function  G,  (for  arbitrary  values  of  U)  the  character  of  the  wave  patterns  changes  at  this 
value  of  i>. 


3.  THEORETICAL  FORMULATION  FOR  SHIP  AND  SEA  WAVES 


To  obtain  the  ship  waves,  the  above  formulation  for  the  elevation  of  a  single  source  oscillating 
at  a  single  frequency  must  be  generalized  to  the  multiple  sources  modeling  the  ship  hull,  and  the  mul¬ 
tiple  frequencies  and  headings  modeling  a  random  sea.  The  first  section  describes  the  procedure  used 
to  model  the  random  sea.  the  following  section  summarizes  the  previously  used  method  [3]  to  gen¬ 
erate  ship  waves  at  zero  forward  speed,  the  third  section  describes  the  modifications  implemented  to 
account  for  the  present  case  of  low  speed,  and  the  last  section  briefly  describes  how  the  wave  slopes 
are  obtained  from  the  elevations. 


The  coordinate  system  used  in  the  present  work  to  calculate  the  ship  and  sea  waves  coincides 
with  that  used  in  the  Ship  Motion  Program  SMP  [14]  and  is  shown  in  Fig.  4.  It  should  be  noted  that 
the  x  axis  is  directed  forward,  from  stern  to  bow.  Thus,  the  ship  forward  speed  in  this  coordinate 
system  corresponds  to  a  flow  in  the  -x  direction  for  the  coordinate  system  for  the  single  pulsating 
source,  shown  in  Fig.  1. 

3.1  Random  Sea  Wave  Elevations 

As  was  done  in  [3],  the  random  sea  elevations  £se(x<  >’*  0  are  approximated  by  a  double  sum 
over  wave  number  km  and  wave  heading  /t„  of  sinusoidal  components,  as  follows 


fS£  =  Re  £  £  -  Re  £  l 


where 


(16) 


m  —  1  n  —  l 


m  =1  #i=l 


rev- 

”*S* 


a,nn  =  V2S(cd„,)Aw  yfG(jin)An  =  ama„ 


S  (o)m )  is  the  value  of  the  frequency  spectrum  at  wave  frequency  u„ 


(17a) 


G(nn)  is  the  value  of  the  spreading  function  at  wave  heading  jxn ,  taken 


1  w* 


relative  to  the  wind  direction  (iw 


:v? 


=  dn  -  a 
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(17b) 
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is  the  wave  heading  measured  with  respect  to  the  x  axis 


Aw  and  An  are  increments  in  os  and  n- 


In  the  present  work,  the  frequency  spectrum  is  taken  to  be  that  of  Pierson-Moskowitz,  while  the 
spreading  function  is  taker,  to  be  the  commonly  used  cosine  squared  function 


IA 


c**. 


2  2 

—  CQS  ft, 
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<  fL  < 


(58) 


As  was  done  in  f3],  five  values  of  o>  are  taken  so  as  to  give  equal  values  of  am,  while  13  equally 
spaced  values  of  ft  between  -r/l  and  tr/2  are  used. 

3.2  Zero  Speed  Ship  Radiation  and  Diffraction  Waves 

Radiation  waves  refer  to  those  generated  by  the  six  oscillatory  modes  of  motion  of  the  ship 
(three  symmetric  or  vertical  modes,  and  three  antisymmetric  or  lateral  modes)  shown  in  Fig.  4,  while 
diffraction  waves  refer  to  the  scattering  of  the  incident  wave  by  the  ship  hull  (one  symmetric,  one 
antisymmetric).  The  resultant  ship  radiation  and  diffraction  waves  £SH  in  the  presence  of  a  random 
sea  are  expressed  in  [3j  by  the  following  triple  sum  over  the  (M  =  5)  sea  frequencies,  ( N  —  13)  sea 
headings,  and  8  ship  oscillation  modes 
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Eft 


ft 


M  N  8  MS  . ,  8 

fw  =  Re  £  £  £  -  Re  £  £  £  7)Jj 

ni  =  1  n  =  1  j  =  1  «i*l  »  =  1  j  =  1 


(19) 


where  T,  (w,  0)  is  the  complex  motion  response  transfer  function  calculated  by  SMP  for  the  six  radia¬ 
tion  modes,  j  =  1,  . . . ,  6,  and  7'7  =  7~g  s  (l.,  0.)  for  the  symmetric  and  antisymmetric  diffraction 
modes,  j  —  7,  8.  The  variable  f :  represents  the  wave  elevation  for  unit  amplitude  for  the  jth  ship 
oscillation  mode  and  is  given  by 


(j  =  Re 


~  ~  Hj  (k ,  6)  (G0  e 
4tt  g  1 


(20) 


where  G0  is  the  Green's  function  for  (J  —  0  given  in  Eq.  (9).  H}  is  the  complex  conjugate  of  the 


Kochin  function  H, 


j 


=  -  j  j  tj  (x’,  y',z')  +w**i  +  r™*\  4S 

s 


(21) 


where  jj  is  the  hull  surface  distribution  of  singularity  strengths  calculated  by  SMP  for  the  jth  radia¬ 
tion  mode,  j  =  1,  ....  6;  x  \  y\  z '  are  respectively  the  x,  >\  z  coordinates  of  the  source  located 
on  the  hull  surface.  The  Kochin  functions  for  the  diffraction  modes  H1  (symmetric)  and  (antisym¬ 
metric)  are  obtained  by  taking  the  singularity  strengths  to  be  a  function  of  wave  heading  0,  as  follows 


77  = 


73  * 


-  ikx  ‘  cosci 


(22a) 


78  =  /  sin  (3  72  e 


ikx '  cos  3 


(22b) 


3.3  Modifications  for  Low  Forward  Speed 


The  above  formulation  must  be  modified  in  a  number  of  ways  in  the  case  of  forward  speed.  In 
view  of  the  approximation  of  low  speed  in  the  derivation  of  G2,  Eq.  (6),  and  the  further  uncertainity 
in  choosing  proper  values  of  the  parameter  p  in  this  equation,  these  modifications  will  be  kept  to  a 
minimum.  The  intent  is  to  model  the  major  features  of  forward  speed,  while  keeping  the  overall  for¬ 
mulation  as  simple  as  possible. 

The  modifications  are  basically  centered  about  the  expression  for  given  in  Eq.  (20).  It  is 
well  known  that  in  the  case  of  forward  speed  U ,  the  radiation  waves  (J  —  1 . 6)  occur  at  the 


encounter  frequency  at  which  is  related  to  wave  frequency  a  and  heading  8  by 


H2e  =  u 


«£/ 

g 


cos  $ 


j  =  1, 


whereas  the  diffraction  waves  (j  =  7,8)  are  not  affected  by  U,  i.e.. 


(23a) 


u€  *  u  j  =  7,  8.  (23b) 

The  second  modification  consists  of  expressing  the  elevation  in  terms  of  G2  instead  of  G0,  by  using 
Eq.  (13)  as  follows 


ML 
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kQ  COS  #0 


G2  (*/') 


(24) 


where  £/'  =  -6'  to  account  for  the  previously  mentioned  difference  in  direction  of  flow  between 
SMP  (Fig.  4)  and  the  single  source  formulation  (Fig.  1). 


Strictly  speaking,  the  Kochin  functions  Hj  and  motion  transfer  coefficients  Tj  for  the  radiation 
modes  j  =  1 . 6  should  be  calculated  for  the  encounter  frequencies  u>e.  Since  there  are  13  head¬ 

ings  and  five  frequencies,  this  would  mean  the  use  of  65  sets  of  Hj  and  Tj.  For  the  sake  of  conveni¬ 
ence,  the  previously  calculated  values  of  Hi  and  T}  for  the  zero  speed  case  [3]  are  also  used  in  the 
present  work.  In  addition  to  the  already  mentioned  approximate  nature  of  the  forward  speed  Green’s 
function  used  here,  the  results  of  [15]  indicate  that  the  motion  coefficients  at  low  speeds  vary  rela¬ 
tively  slowly  with  speed.  Also,  as  shown  later,  the  diffraction  waves  (j  —  7,  8),  for  which  the 
singularity  distribution  for  zero  speed  is  correct,  tend  to  be  larger  than  the  radiation  waves.  Thus,  the 
present  approximation  qualitatively  preserves  the  major  effects  of  forward  speed. 

Combining  Eqs.  (23),  (24),  and  the  preceding  remarks,  the  low  speed  equivalents  for  the  zero 
speed  single  mode  and  total  ship  waves,  respectively  given  by  Eqs.  (20)  and  (19),  become 
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m  —  1  n  —  1  J  =  1 


i  —  i  n  —  1 


j  - 1 


3.4  Wave  Slopes 

As  shown  in  Eq.  (14),  wave  slopes  are  conveniently  obtained  by  multiplying  the  elevations  by 
the  wavenumbers  in  the  desired  horizontal  direction.  The  wave  slopes  in  the  x  and  y  directions  for 
the  sea  waves  f5£  (Eq.  16),  ship  waves  (Eq.  26).  and  superposed  waves  are  simply  given  by 
the  following  sums 
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where 
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cmj 


g  +2U'  COS 


Cy  »s  given  by  g  sin(7  ~  d^j)  +  2w„,y  U'  sin  7  =  0 


is  defined  by  Eqs.  (23a)  and  (23b). 


4.  COMPUTER  CALCULATIONS 
4.1  Computer  Program 

The  previously  used  code  for  zero  speed  (3]  has  been  revised  to  account  for  low'  forward  speed. 
The  program  accepts  from  one  input  file  data  on  the  speed  and  direction  of  the  wind,  the  frequency 
spectrum  to  be  used,  the  number  of  frequency  components,  the  types  of  waves  to  be  calculated  (ship, 
sea,  or  total  of  ship  and  sea;  elevation,  x  slope,  or  y  slope)  and  the  ship  speed.  In  this  same  file  are 
specified  data  for  the  rectangular  grid  over  which  the  waves  are  calculated:  its  origin,  and  the 
number  and  spacing  of  the  points  along  the  (a  ’,  y')  axes,  which  are  inclined  at  angle  a  to  the  (a,  y) 
system  (see  Fig.  5).  As  a  second  input  file,  the  program  accepts  the  Kochin  functions  Wy(u >.  0.  j3) 
and  transfer  functions  r,(cj.  (3)  calculated  by  SMP  (14). 

The  program  then  performs  the  calculations  for  the  wave  elevations  and  slopes  given  by  Eqs. 
(16)-(29).  The  calculated  waves  are  saved  as  an  output  plot  file  from  which  three-dimensional  or 
contour  plots  may  be  generated.  Examples  of  these  plots  are  given  in  the  following  chapter. 
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4.2  Adaptation  of  Formulation  to  SMP  Results 

SMP  calculates  the  functions  Hj  and  7}  for  30  discrete  values  of  «  between  0.2  and  2.4  rad/sec, 
and  13  equally  spaced  values  of  6  and  £  between  0  and  180  degrees.  Thus,  it  is  of  interest  to 
describe  the  procedure  used  to  model  the  general  case  of  values  of  w,  j3,  and  $  which  may  not  coin¬ 
cide  with,  or  fail  outside  the  range  of,  the  SMP  results.  This  procedure  is  described  in  detail  in  [3], 
and  is  briefly  summarized  here. 

In  order  to  more  accurately  model  the  lower  sea  states,  where  the  higher  w’s  are  of  interest,  the 
two  highest  frequencies  calculated  by  SMP  were  changed  from  2.2  and  2.4  to  2.4  and  2.8  rad' sec, 
respectively.  This  essentially  consisted  of  changing  one  DATA  statement  in  SMP.  The  values  of  Hj 
and  T,  are  approximated  by  those  corresponding  to  the  SMP  values  of  w  which  come  closest  to  the 
actual  values. 

The  wave  headings  $  are  taken  to  be  the  13  equally  spaced  headings  between  0  and  180  degrees 
used  in  SMP.  Noting  that  the  energy  of  the  wind  is  spread  90  degrees  on  either  side  of  the  wind 
direction  {see  Eq.  {18)),  the  SMP  formulation  is  most  conveniently  used  when  the  wind  direction  fiw 
is  90  degrees,  i.e.,  in  this  case  the  wave  headings  of  interest  are  simply  the  0  to  180  degree  cases 
considered  by  SMP.  For  other  wind  directions,  not  all  13  values  considered  by  SMP  are  used,  and 
negative  values  of  $  are  involved.  In  these  cases,  use  is  made  of  the  fact  that  for  the  vertical  modes 
(j  -  1,3.  5,  7)  Tj(u.  -0)  -  T.(u>,  0)  while  for  the  lateral  modes  (j  -  2, 4,  6,  8) 

Tj(u.  -&)  m  -  Tj{u,  &),  i.e., 

7}<«,  -J3)  =  <-iy  +  1  Tj(u>,  0)  j  =  1,2,...,  7.  8.  (30) 

The  Koch  in  functions  H}  are  defined  only  for  13  positive  values  of  6  evenly  spaced  between  0 
and  180  degrees.  For  arbitrary  positive  values  of  6,  the  value  of  Hj  is  obtained  by  simple  interpola¬ 
tion.  For  negative  values  of  0,  Hj  is  obtained  by  noting  that  the  ship  waves  are  symmetrical  with 
respect  to  0  for  the  vertical  modes  and  antisymmetric  for  the  lateral  modes,  i.e., 

Hj(-6)  =  (-iy  +  !  Hj(6)  j  =  I,  2 . 7,  8.  (31) 


4.3  Computer  Runs 

In  view  of  the  approximations  made  in  the  present  work,  it  is  not  felt  that  the  theory  can  effec¬ 
tively  differentiate  ship  hulls  of  different  shape.  Thus,  all  runs  were  made  using  the  DE1006  as  the 
sample  hull.  This  ship  has  a  length  L  of  308  ft,  a  beam  B  of  35.9  ft,  and  a  draft  H  of  12.1  ft. 

Three  Sea  States  were  considered:  2,5,  and  8,  which  respectively  correspond  to  wind  speeds  of 
12,  24,  and  42  knots.  Table  2  shows  the  five  frequencies  oj  and  associated  wavelengths  X  chosen  to 
model  each  Sea  State.  This  table  shows  that  the  ranges  of  X’s  for  Sea  State  2,  5,  and  8  are  respec¬ 
tively  smaller  than,  include,  and  larger  than  the  ship  length  L.  Most  of  the  runs  were  made  for  a 
wind  heading  of  45  degrees,  which  corresponds  to  quartering  seas.  To  illustrate  the  effect  of 
heading  on  the  ship  waves,  runs  were  also  made  where  heading  varied  between  0  and  180  degrees. 
Two  values  of  (/  were  considered  in  the  computer  runs:  U  —  0  and  U  —  10  knots  (16.88  ft/s), 
which  gives  the  Froude  number  F„  =0.17  where 


(32) 
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5. 


CALCULATED  SEA  AND  SHIP  WAVE  PATTERNS 


5.1  Random  Sea  Elevations  and  Slopes 

Figures  6  and  7  show  plots  of  the  random  sea  elevations  and  slopes.  These  are  given  for  illus¬ 
trative  purposes  as  well  as  for  comparison  with  the  ship  generated  waves  shown  later  in  this  chapter. 
Figures  6a,  b,  and  c  respectively  show  three-dimensional  plots  of  the  elevations  for  Sea  States  2,  5, 
and  8,  while  Figs.  6d-f  show  the  corresponding  contour  plots  for  these  Sea  States.  Figures  7a  and  b 
respectively  show  three-dimensional  plots  of  the  x  and  y  slopes  for  Sea  State  5,  while  Figs.  7c  and  d 
show  the  corresponding  contour  plots. 

The  w'ind  direction  (3W  is  taken  to  be  45  degrees,  the  quartering  sea  case.  The  wave  patterns 
are  calculated  for  a  square  x  -y  grid  given  by 

0  <  .t .  y  <  4  Xo  (33) 

where  Xq  is  the  longest  wavelength,  corresponding  to  the  fundamental  frequency  oj0,  shown  for  each 
Sea  State  in  Table  2. 

5.2  Effect  of  vl;  on  Ship  Waves 

It  was  previously  shown  in  Fig.  3  the  manner  in  which  the  wave  elevations  for  a  single  source 
become  singular  as  the  upper  limit  of  v ,  vu,  approaches  0.5.  As  a  result,  it  was  pointed  out  that  vv 
is  taken  equal  to  0.25  in  most  of  the  ship  wave  calculations.  Figures  8a  and  b  show,  for  the  case  of 
Sea  State  8,  j 3W  —  45  degrees,  f/  =  10  knots,  three-dimensional  plots  of  the  ship  waves  fSH  and  the 
superposed  ship  and  sea  waves  for  the  square  grid  defined  by  Eq.  (33)  for  a  choice  of  vv  -  0. 
Figs.  8c-f  show  corresponding  pairs  of  plots  for  choices  of  vv  =  0.25  and  0.45. 

Figures  8b  and  d  show  that  for  vv  <  0.25,  there  is  little  visual  difference  between  the  super¬ 
posed  waves  tr  and  the  sea  waves  £SE  shown  in  Fig.  6c  for  Sea  State  8.  When  vv  increases  to 
0.45,  Figs.  8e  and  f  show  that  now  show  features  which  are  clearly  those  of  Results  for 
vv  =  0.49,  which  are  not  presented,  show  that  f SH  becomes  so  iarge  that  it  is  virtually  indistinguish¬ 
able  from  Tf.  Thus,  similar  to  the  single  source  case,  the  ship  waves  can  be  made  arbitrarily  large 
by  choosing  vv  arbitrarily  close  to  0.50.  In  all  subsequent  ship  waves  shown  in  this  chapter,  the 
intermediate  value  of  0.25  is  used. 

5.3  Effect  of  Forward  Speed  U 

Figures  9a-d  show  three-dimensional  plots  of  the  ship  wave  elevations  $SH  for  Sea  State  5, 
=  45  deg,  as  U  increases  from  0  to  20  knots  (33.76  ft/s).  The  elevations  are  plotted  for  a  rec¬ 
tangular  grid  which  now  shows  the  waves  aft  as  well  as  forward  of  the  ship 

-  4  Xo  <  x  <  4  Xq,  0  <  y  <  4  Xq.  (34) 

The  figure  shows  how  the  upstream  wpve-s  tend  to  cluster  and  the  downstream  waves  tend  to  disperse 
beyond  the  end  of  the  grid  as  U  increases.  Due  to  the  adopting  of  the  upper  limit  of  0.25  on 
v  -  U&ig ,  Figs.  9c  and  d  show  that  there  is  relatively  little  change  in  both  the  pattern  and  amplitude 
of  the  ship  waves  for  J  2:  6  knots  (10.13  ft/s). 


5.4  Effect  of  Wind  Heading  /3VV 


Figures  10a  and  b  show  three-dimensional  plots  of  the  ship  waves  for  Sea  State  5,  U  -  0,  over 
the  rectangular  grid  defined  in  Eq.  (34),  for  values  of  equal  to  45  and  135  deg.  Figures  11a  and 
b  show  corresponding  results  for  U  =  10  knots.  Figure  10  shows  that  for  U  =  0  the  strongest 
waves  are  those  which  are  diffracted  approximately  at  the  same  angle  as  (3W ,  and  there  is  a  certain 
symmetry  between  the  plots  for  and  1 80-/3^.  Due  to  the  upstream  clustering  and  downstream 
dispersion  of  the  waves,  Fig.  1 1  shows  that  no  such  symmetry  exists  in  the  forward  speed  case. 
Results  for  (3W  =  0  and  180  deg  show  a  similar  trend.  Thus,  for  wind  headings  of  0  and  45  deg, 
most  of  the  waves  are  clustered  upstream.  On  the  other  hand,  for  the  corresponding  wind  headings  of 
180  and  135  deg,  the  ship  waves  fill  the  entire  calculation  grid. 

5.5  Decomposition  of  Ship  Waves  Into  Radiation  and  Diffraction  Components 

For  Sea  State  5  with  (3W  —  45  deg,  U  =  10  knots,  and  the  square  calculation  grid  given  by  Eq. 
(33),  Figs.  12a  to  c  respectively  show  three-dimensional  plots  of  the  resultant  ship  waves  and 
the  radiation  and  diffraction  components,  £SHR  and  £shd .  Recall  that  £SHR  refers  to  the  sum  over 
j  =  1,  . . .  ,  6  and  f SHD  refers  to  the  sum  over  j  =  7,  8  of  the  expression  for  given  by  Eqs.  (19) 
or  (26). 

This  figure  clearly  shows  that  the  major  contribution  to  the  resultant  ship  waves  is  made  by  the 
diffraction  component.  Similar  results  hold  for  other  Sea  States  and  values  of  U.  This  tends  to  jus¬ 
tify  the  simplification  discussed  in  Section  3.3  of  using  the  zero  speed  Kochin  functions,  which  is 
appropriate  only  for  the  diffraction  waves.  This  figure  also  shows  the  expected  trend  of  the  diffrac¬ 
tion  waves  being  more  directional  than  the  radiation  waves. 

5.6  Comparison  of  Superposed  and  Sea  Waves 

For  =  45  deg,  U  =  0,  and  the  square  calculation  grid  given  by  Eq.  (33),  Fig.  13a  shows 
the  three-dimensional  plot  of  the  superposed  wave  elevations  for  Sea  State  2,  while  Figs.  13b  and 
c  show  the  corresponding  plots  for  t/  =  10  knots  for  Sea  States  2  and  5.  Figures  14a  to  c  show  con¬ 
tour  plots  of  these  same  elevations.  Figures  15a  and  b  respectively  show  three-dimensional  plots  of 
the  x  and  y  slopes  of  the  superposed  waves  for  Sea  State  5  for  U  =  0,  while  Figs.  15c  and  d  show 
the  plots  for  the  U  =  10  knots  case.  Finally,  Figs.  16a  and  b  show  the  corresponding  contour  plots 
of  the  U  =  0  case  given  in  Figs.  15a  and  b. 

The  superposed  waves  in  the  above  figures  may  be  compared  to  the  ocean  waves  shown  in  Figs. 
6  and  7.  In  the  case  of  the  three-dimensional  plot  for  (7=0,  Fig.  13a  shows  that  even  in  the  case  of 
Sea  State  2,  for  which  diffraction  effects  dominate,  the  superposed  waves  are  only  moderately  dif¬ 
ferent  from  the  corresponding  sea  waves  shown  in  Fig.  6a.  The  plots  of  the  superposed  waves  for 
Sea  States  5  and  8  are  virtually  indistinguishable  from  the  corresponding  sea  waves,  and  hence  have 
not  been  shown.  In  the  case  of  U  —  10  knots,  for  which  the  ship  waves  are  larger,  Fig.  13b  shows 
that  in  the  case  of  Sea  State  2,  the  superposed  waves  are  completely  different  from  the  sea  waves 
while  the  differences  arc  small  but  now  perceptible  for  Sea  State  5.  The  same  observations  hold  in 
comparing  the  contour  plots  of  the  superposed  waves  in  Figs.  14a  to  c  with  the  corresponding  sea 
waves  shown  in  Figs.  6d  and  c.  In  all  cases,  the  trend  for  Sea  State  8  is  similar  to  that  of  Sea  State 
5. 


In  view  of  the  small  differences  between  the  superposed  and  sea  wave  elevations  for  the  higher 
Sea  States,  it  is  of  interest  to  compare  the  wave  slopes.  In  the  case  of  three-dimensional  plots  for  Sea 
State  5.  Figures  15a  and  b  show  that  the  x  and  y  slopes  of  the  superposed  waves  are  only  slightly  dif¬ 
ferent  from  the  corresponding  sea  wave  slopes  shown  in  Figs.  7a  and  b.  On  the  other  hand,  for 
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U  —  10  knots,  Figs.  15c  and  d  show  that  there  are  substantial  differences  in  the  case  of  the  x  slope 
and  moderate  differences  in  the  case  of  the  y  slope.  The  larger  differences  for  the  x  slope  are  due  to 
the  previously  mentioned  wave  clustering  in  the  upstream  direction,  which  leads  to  shorter 
wavelengths.  The  near  lack  of  difference  between  the  three-dimensional  plots  of  the  wave  slopes  for 
U  —  0  makes  it  of  interest  to  further  consider  the  corresponding  contour  plots.  Figure  i6a  shows 
that  the  differences  are  very  small  for  the  x  slope,  but  noticeable  for  the  y  slope.  While  not  shown,  a 
similar  difference  exists  in  the  y  slope  at  U  —  0  for  Sea  State  8. 

6,  SUMMARY 

A  previously  derived  formulation  for  calculating  ship  generated  radiation  and  diffraction  waves 
has  been  generalized  to  the  case  of  low  forward  speed.  One  modification  consists  of  replacing  the 
previously  used  far  field  Green’s  function  for  zero  speed  by  a  low  forward  speed  function,  for  which 
all  terms  involving  the  square  of  speed  are  considered  negligible.  It  is  shown  that  this  new  function 
retains  the  previous  simplicity  of  one  set  of  waves  everywhere,  with  the  difference  that  the  wave  pro¬ 
pagation  direction  no  ionger  coincides  with  the  spatial  direction  of  the  field  point.  The  manner  in 
which  the  wavelengths  decrease  and  wave  amplitudes  increase  is  shown  as  the  singular  limit  of  0.5  is 
approached  for  the  parameter  v  -  U w/g ,  where  U  is  the  forward  speed  and  o>  is  the  circular  fre¬ 
quency.  For  an  upper  bound  on  v,  vv,  equal  to  0.25,  the  maximum  wave  amplitude  is  approximately 
three  times  as  large  as  that  of  the  zero  speed  case.  A  second  modification  to  the  zero  speed  case  con¬ 
sists  of  replacing  the  sea  wave  frequencies  by  frequencies  of  encounter  for  the  radiation  modes. 

A  series  of  three-dimensional  and  contour  plots  are  given  for  the  sea,  ship,  and  superposed  sea 
and  ship  wave  elevations  and  slopes  for  the  DE1006  ship  in  the  presence  of  Sea  States  2,  5,  and  8.  It 
is  again  shown  that  by  choosing  vv  arbitrarily  dose  to  0.5  the  ship  waves  can  always  be  made  to 
dominate  the  sea  waves.  The  tendency  for  the  upstream  waves  to  cluster  and  the  downstream  waves 
to  disperse  as  U  increases  is  illustrated  by  a  series  of  three-dimensional  plots.  For  U  ~  0,  there  is  a 
certain  near  symmetry  between  the  ship  waves  for  wind  headings  @w  and  180-0#/.  For  U  -  10 
knots,  no  such  symmetry  exists  due  to  the  above  upstream  clustering  and  downstream  dispersion.  For 
both  <7  =  0  and  10  knots,  the  ship  diffraction  waves  tend  to  dominate  the  radiation  waves. 

In  comparing  the  superposed  ship  and  sea  waves  with  the  corresponding  sea  waves  alone,  the 
plots  for  zero  forward  speed  show  little  difference  in  the  elevations  and  x  slopes  of  these  waves  for 
Sea  States  5  and  8,  and  moderate  differences  in  the  case  of  Sea  State  2.  It  is  only  in  the  contour  plots 
of  the  y  slope  that  appreciable  differences  can  be  detected  between  the  superposed  and  sea  waves  at 
the  higher  Sea  States^  In  the  case  of  V  =  10  knots,  for  which  the  ship  waves  are  larger,  there  now 
exist  small  (but  perceptible)  differences  for  the  wave  elevations,  large  differences  in  the  x  slope,  and 
moderate  differences  in  the  y  slope  at  the  higher  Sea  States. 
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Fig.  1  —  Definition  of  Grekas’s  coordinate  system  for  single  source 


16 


Fig.  2  (Continued)  —  Sample  wave  patterns  for  single  source  for  different  cases  of  U 


Fig.  3d  —  vy  —  0.40 
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Fig.  3e  —  vv  =  0.49 


Fig.  3  (Continued)  —  Sample  wave  patterns  for  single  source  for  V  —  30  ft/s,  different  values  of 
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Fig.  6d  —  Contour,  Sea  State  2 

Fig.  6  (Continued)  —  Plots  of  random  sea  wave  elevations,  =45  degrees 


Fig.  8a  —  Ship  waves,  vv  -  0 
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Fig,  8b  —  Superposed  waves,  vv  —  0 


Three-dimensional  plots  of  ship  and  superposed  wave  elevations  for  Sea  State  8. 
-  45  degrees,  different  values  of  vy 


Fig.  8d  —  Superposed  waves,  vv  =  0.25 
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n. snued)  —  Three-dimensional  plots  of  ship  and  superposed  wave  elevations  for  Sea  State  8, 
Pw  =  45  degrees,  different  values  of  vjj 


Fig.  9  (Continued)  —  Three-dimensional  plots  of  ship  wave  elevations  for  Sea  State  5, 

—  45  degrees,  different  values  of  U 


Fig.  10  —  Three-dimensional  plots  of  ship  wave  elevations  for  Sea  State  5, 
V  —  0,  different  values  of  0K- 
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Fig.  1  la  —  0W  —  45  degrees 
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Fig.  lib  —  —  335  degrees 


Fig.  11  —  Three-dimensional  plots  of  ship  wave  elevations  for  Sea  State  5, 
U  =  10  knots,  different  values  of 


Three-dimensional  plots  of  superposed  wave  elevations,  =  45  degrees 


Fin.  !3c  —  Sea  State  5,  U  —  10  knots 


13  (Continued)  —  Three-dimensional  plots  of  superposed  wave  elevations.  / =  45  degrees 
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Fig.  14b  —  Sea  State  2,  U  =  iO  knots 
—  Contour  plots  of  superposed  wave  elevations,  Pw  ~  45  degrees 


Fig.  14c  —  Sea  State  5,  U  =  10  knots 

Fig.  14  (Continued)  —  Contour  plots  of  superposed  wave  elevations,  0W  =  45  degrees 


Fig.  15  (Continued)  —  Three-dimensional  plots  of  superposed  wave  slopes 
for  Sea  State  5,  (3W  =45  degrees 


Fig.  15d  —  }> -slope,  U  —  10  knots 


Fig.  15  (Continued)  —  Three-dimensional  plots  of  superposed  wave  slopes 
for  Sea  State  5,  0W  =45  degrees 


41 


Fig.  16b  —  y -slope 

Fig.  16  —  Contour  plots  of  superposed  wave  slopes  for  Sea  State  5,  U  -  0,  Qw  -  45 


42 


